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1 Introduction 


Knowledge of oceanic irradiance is critical for realistic simulation of ocean biogeochem- 
istry. Since irradiance affects ocean thermodynamics and physiological processes of 
biota in the oceans, it may affect the abundances and distributions of phytoplankton, 
and potentially affect the uptake of carbon dioxide. 

In order to understand the effects of irradiance on phytoplankton abundances and pri- 
mary production in the global oceans, we require a reasonably realistic simulation of 
the irradiance reaching the ocean surface, and propagating through the water column. 
The simulation must explicitly include the effects of atmospheric optical constituents 
on surface irradiance, processes occurring at the air-sea interface, and the results of 
these effects in the water column taking into account the optical properties of the water, 
phytoplankton, and dissolved detrital materials. In this paper we describe atmospheric 
and oceanic radiative transfer models, combined as the Ocean-Atmosphere Radiative 
Model (OARM) that provide a reasonably complex and realistic representation of at- 
mospheric radiative processes and the irradiance availability in the oceans. 


2 Ocean- Atmosphere Radiative Model (OARM) 

2.1 Atmospheric Component 

The atmospheric radiative model is based on the Gregg and Carder (1990) spectral 
model for clear skies, and relies on Slingo (1989) for spectral cloud transmittance. The 
clear sky model of Gregg and Carder (1990) was derived from Bird and Riordan (1986), 
but was limited to the spectral range of photosynthetically available radiation (PAR), 
defined as 350-700 nm. It also contained increased spectral resolution and marine 
atmospheric and surface reflectance conditions. It compared within ±6.6% root mean 
square (RMS) with surface observations of spectral irradiance and ±5.1% RMS with 
integrated PAR (Gregg and Carder, 1990). This model is extended for OARM from the 
PAR spectral domain to the entire solar spectrum, from 200 nm to 4 um, representing 
>99% of the total solar irradiance impinging on the top of the atmosphere. For compu- 
tational efficiency, the spectral resolution is degraded from 1 nm used in Gregg and 
Carder (1990) to a variable resolution appropriate for the spectral absorbing properties 
of the major atmospheric optically active gases, specifically ozone, water vapor, oxy- 
gen, and carbon dioxide (Table 1). The spectral resolution is fixed at 25 nm for the 
PAR range, which is the region of interest for phytoplankton photosynthesis. 

A modification to the clear sky model is the use of spectral surface reflectance in the 
presence of sea foam based on observations by Frouin et aL (1996), and confirmed by 
M oore et al. (1998). This modification accounts for the fact that some of the reflec- 
tance by sea foam is due to bubbles located slightly below the surface. This produces 
some absorption by sea water before reflecting the light out of the water. Since longer 
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Table 1. Atmospheric optical data for OARM. 1 is wavelength (nm), F 0 is extraterrestrial 
irradiance (W/m2), t r is Rayleigh optical thickness (dimensionless), a oz is ozone 
absorption coefficient (cm 1 ), a wv is water vapor absorption coefficient (cm 1 ), a 0 2 is 
oxygen absorption coefficient (cm' 1 ), and a C0 2 is carbon dioxide absorption coefficient 
(cm -1 ). A data file is available online at salmo.gsfc.nasa.gov, 
pub/outgoing/oarm/atmo25.dat. 


1 

F 0 

tr 

3 oz 

&WV 


& eo 2 

250 

16.5280 

2.8229 

80.8836 

0.0000 

0.0000 

0.0000 

325 

41.0035 

0.8742 

0.3584 

0.0000 

0.0000 

0.0000 

350 

12.7915 

0.5925 

0.0016 

0.0000 

0.0000 

0.0000 

375 

27.1180 

0,4768 

0.0001 

0.0000 

0.0000 

0.0000 

400 

35.1935 

0.3643 

0.0000 

0.0000 

0.0000 

0.0000 

425 

42.1795 

0.2833 

0.0008 

0.0000 

0.0000 

0.0000 

450 

49.1940 

0.2238 

0.0040 

0.0000 

0.0000 

0.0000 

475 

49.7700 

0.1791 

0.0121 

0.0000 

0.0000 

0.0000 

500 

48.3945 

0.1451 

0.0279 

0.0000 

0.0000 

0.0000 

525 

46.4935 

0.1189 

0.0529 

0.0000 

0.0000 

0.0000 

550 

46.5160 

0.0983 

0.0832 

0.0000 

0.0000 

0.0000 

575 

46.1380 

0.0820 

0.1140 

0.0148 

0.0000 

0.0000 

600 

44.1230 

0.0690 

0.1176 

0.0851 

0.0000 

0.0000 

625 

42.0468 

0.0584 

0.0967 

0.0016 

0.0058 

0.0000 

650 

39.5245 

0.0498 

0.0643 

0.0408 

0.0000 

0.0000 

675 

37.8730 

0.0427 

0.0416 

0.0009 

0.2333 

0.0000 

700 

18.7465 

0.0382 

0.0261 

0.1074 

0.1299 

0.0000 

725 

66.8775 

0.0319 

0.0132 

0.1462 

0.0000 

0.0000 

775 

59.9900 

0.0244 

0.0072 

0.0065 

1.1176 

0.0000 

850 

102.2358 

0.0168 

0.0020 

0.0568 

0.0001 

0.0000 

950 

84.7899 

0.0108 

0.0000 

4.9142 

0.0000 

0.0000 

1050 

67.8772 

0.0072 

0.0000 

0.0279 

0.0002 

0.0002 

1150 

56.0762 

0.0050 

0.0000 

17.5902 

0.0000 

0.0000 

1250 

46.1106 

0.0032 

0.0000 

0.0479 

0.0028 

0.0005 

1350 

37.3780 

0.0026 

0.0000 

115.3197 

0.0000 

0.0003 

1450 

31.8313 

0.0020 

0.0000 

95.2031 

0.0000 

0.0023 

1550 

27.2524 

0.0015 

0.0000 

0.0364 

0.0000 

0.0014 

1650 

22.7691 

0.0012 

0.0000 

0.0082 

0.0000 

0.0010 

1750 

18.6144 

0.0009 

0.0000 

2.0008 

0.0000 

0.0001 

1900 

27.9521 

0.0007 

0.0000 

228.9211 

0.0000 

0.0087 

2200 

33.6300 

0.0004 

0,0000 

0.2550 

0.0000 

0.0118 

2900 

32.5775 

0.0001 

0.0000 

314.7111 

0.0000 

0.0256 

3700 

7.2635 

0.0000 

0.0000 

0.9158 

0.0000 

0.0003 
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wavelengths are much more strongly absorbing than shorter ones, the net effect is re- 
duction of foam reflectance at longer wavelengths. Based on Frouin’s data (Frouin et 
al., 1996) foam reflectance is treated as 

F(l) = a + a ln[T (1)] + a ln[T (l)] 2 + a InfT (l)] 3 ; 1< 900nm (1) 

0 1 w 2 w 3 w 

where F is a factor to adjust foam reflectance for spectral dependence, 1 is wavelength 
in nun, T w is the transmittance of water 

T = exp [a (1) + 0,5b (l)J (2) 

w w w 

where a w and b w are the spectral absorption and scattering coefficients of seawater, 
respectively. The 0.5 factor for b w represents the backscattered portion. In Eq, 1, aO, 
al, a2, and a3 are empirical coefficients to produce a best fit to the published data of 
Frouin et al. (1996), and are 0.9976, 0.2194, 0.0554, and 0.0067, respectively. 

At and above 900 nm, a w is nearly completely absorbing, and a new parameterization is 
used 


F(l) = b 0 + b ; l + b,l 2 + b 3 l 3 ; 1 > 900 nm (3) 

where b 0 , b 15 b 2 , and b 3 are defined similarly to their counterparts in Eq. 1, and are 
5.026, -0.0114, 9.552xl0- 6 , and 2.698x10". respectively. 

Total foam reflectance r f now includes spectral dependence through 

r.tW.h = r f (W)F(l) (4) 

where W is the wind speed (m s _1 ) and r f (W) is defined as in Gregg and Carder (1990). 
Total surface reflectance is the sum of the specular and foam contributions, which now 
have a spectral dependence 


r (W,l) = r (W) + r (W,l) (5) 

d dsp f 

r (W,l) = r (W) + r (W,l) (6) 

s ssp f 

where r d and r s represent the direct and diffuse surface reflectance components, and the 
subscript sp represents specular reflectance (Gregg and Carder, 1990). 

Under cloudy skies, OARM utilizes the cloudy sky model of Slingo (1989). This model 
computes irradiance transmittance through clouds based on the liquid water path and 
the cloud droplet effective radius using a Delta-Eddington approximation of the two- 
stream approach. It is used for shortwave cloud forcing in the National Center for 
Atmospheric Research Community Climate Model 3 (Kiehl et al., 1998). 
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In our implementation of the S lingo (1989) model, gaseous absorption by ozone, oxy- 
gen, carbon dioxide, and water vapor, occurs before cloud transmittance, and Earth 
curvature effects (Kasten, 1966) are accounted for. The cloud droplet radius is set to a 
constant 1 1.8 um, which represents the mean oceanic value observed in a global survey 
(Han et al„ 1994). 

The clear and cloudy sky models track two irradiance streams to the ocean surface: 
direct and diffuse (Fig. 1). The magnitude, spectral, and directional effects of clouds 
on irradiance reaching the ocean surface are illustrated in Fig. 2 for a liquid water path 
of 80 g nr 2 . This water path value produces about one-half the irradiance of clear skies 
for the same solar zenith angle (30°), atmospheric optical gases and molecular and par- 
ticulate scattering. Clouds clearly change the amount of irradiance reaching the ocean 
surface, but also directional quality by apportioning a larger percentage to diffuse rather 
than direct, which predominates under clear skies. The spectral effects of clouds are 
shown more clearly in Table 2, where proportions of different parts of the visible spec- 
trum are listed with and without clouds. Cloudy sky blue spectra (350-400 nm) repre- 
sent a larger proportion of total PAR than clear sky spectra (20.5%), while cloudy sky 
red spectra (650-700 nm) represent a smaller proportion (-9.4%). These results con- 
firm observations by Siegel et al. (1999). 



Figure 1 . Depiction of the pathways of irradiance under clear and cloudy skies, and 
in the oceans. The sizes of the arrows indicate the relative proportions of direct (E d ) 
and diffuse (E s ) irradiance for clear skies and cloudy skies. Some of the surface 
irradiance is reflected off the sea surface (1-r). These pathways continue into the 
ocean where an additional diffuse upwelling (E n ) path exists. 
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Figure 2. Spectral surface irradiance just below the sea surface (after spectral surface re- 
flectance) for clear skies and cloudy skies. The cloudy sky simu lation represents the effects 
of 80 g m 2 liquid water path, which produces about half the total surface solar irradiance as 
the clear sky model for the same solar zenith angle and atmospheric optical properties. 


Table 2. Comparison of proportion of blue irradiance to 
red irradiance under clear and cloudy skies. Cloudy sky 

used a liquid 

water path of 80 g mC 

l, which produces a 

reduction just under the sea surface of about half that for 

clear skies. 




Percent of PAR 


350-400nm 

650-700 nm 

Clear 

11.2% 

17.0% 

Cloudy 

13.5% 

15.4% 

Difference 

20.5% 

-9.4% 


(3 ARM requires external cloud properties (cloud cover and liquid water path), as well 
as surface pressure, wind speeds, relative humidity, precipitable water, and ozone. 
Climatologies of cloud variables were created using a 5-year record of International 
Satellite Cloud Climatology Project (ISCCP)-D2 monthly mean data made available by 
the NASA/Goddard Earth Sciences Distributed Active Archive Center (GES-DAAC). 
Climatologies of surface pressure, wind speeds, relative humidity, precipitable water 
were created using National Center for Environmental Prediction (NCEP) Reanalyses 
for the years 1978-2000 (Kalnay et al„ 1996). Climatological ozone was derived from 
a 20.5 year record beginning in 1978 of several Total Ozone Mapping Spectrometer 
(TOMS) sensors, again made available by the GES-DAAC. 
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2.2 Oceanic Component 


The oceanic portion of OARM tracks three irradiance streams in the water column: 
downwelling direct, downwelling diffuse, and up welling diffuse (Fig. 1). A complete 
mathematical description of the oceanic component is described in the Appendix. How- 
ever. the highlights of the model are described in the following. 

The direct downwelling stream contains the irradiance directly transmitted by the sun, 
traversing the air-sea interface, and proceeding forward at an angle described by the 
solar zenith angle modified by the refractive index of seawater. Each scattering and 
absorbing event in the water column removes irradiance from the direct downwelling 
stream. 

Whereas the downwelling direct irradiance receives no contributions in the water col- 
umn and steadily decreases as the result of absorption and scattering processes, the 
downwelling diffuse irradiance gains forward-scattered downwelling direct irradiance 
and backscattered upwelling diffuse irradiance in the water column. This irradiance 
stream travels along a path defined as the average cosine for downwelling diffuse irra- 
diance, This is considered constant in OARM at 0.83 (Sathyendranath and Platt, 1988). 

The diffuse upwelling stream is composed entirely of backscattered direct and diffuse 
downwelling irradiance. It is derived from the two-stream model of Aas (1987), with 
contributions from both the downwelling direct and downwelling diffuse explicitly in- 
cluded, similarly to Ackleson et al. (1994). The diffuse upwelling stream can represent 
10% of the total at some wavelengths. We track only one scattering event for each of 
the oceanic layers for computational convenience. This approximation is effective for 
the open ocean domain, but would become erroneous in waters containing very large 
amounts of optically active constituents. 

Application of the oceanic portion of OARM requires that the absorption and scattering 
properties of the water column be specified. In this effort we include fixed contribu- 
tions by seawater and chromophoric dissolved organic matter (CDOM), and variable 
contributions by five different phytoplankton groups. Water absorption and scattering 
properties are taken from Smith and Baker (1981) for the range 200-800 nm, Circio and 
Petty (1951) for the range 800nm-2.5mm, and Maul (1985) for the range 2. 5-4. Omm. 
A reanalysis of seawater absorption for 380-729 nm (Pope and Fry, 1997), is substi- 
tuted for the Smith and Baker (1981) estimates (Table 3). CDOM absorption is taken 
from Carder et a I, (1989), which includes a humic acid component and a fulvic acid 
component. 

Phytoplankton specific absorption coefficients are derived by taking spectra from the 
published literature and normalizing to the absorption at 440 nm [a* p (440)] to obtain 
normalized specific absorption spectra [a* p (l)] N for each of the five phytoplankton groups 
(see Appendix). The resulting spectra exhibit substantial inter-group variability (Fig. 
3). 
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Table 3. Spectral seawater absorption (av) and total scattering (b w ) coefficients (m 1 ). 
Derived from Smith and Baker (1981) for the spectral ranges 200-370 nm and 730-800 
nm, Pope and Fry (1997) for the range 380-720 nm, Circio and Petty (1951) for the 
range 800nm-2.5nm, and Maul (1985) for the range 2.5-4nm. 


I (nm) 

a w 

bw 

250 

0.6112 

0.0567 

325 

0.0762 

0.0187 

350 

0.0461 

0.0135 

375 

0.0182 

0.0100 

400 

0,0063 

0.0076 

425 

0.0051 

0.0058 

450 

0.0083 

0.0045 

475 

0.0119 

0.0036 

500 

0.0215 

0.0029 

525 

0.0407 

0.0023 

550 

0.0550 

0.0019 

575 

0.0849 

0.0016 

600 

0.1995 

0.0014 

625 

0.2850 

0.0012 

650 

0,3512 

0.0009 

6 75 

0,4559 

0.0007 

700 

0.6433 

0.0007 

725 

j .4449 

0.0006 

775 

2.3900 

0.0004 

850 

3.7382 

0.0002 

950 

27.4805 

0.0000 

1050 

19.3470 

0.0000 

1150 

67.1800 

0.0000 

1250 

94.9976 

0.0000 

1350 

363.1256 

0.0000 

1450 

1118.6070 

0.0000 

1550 

944.8757 

0.0000 

1650 

519.5995 

0.0000 

1750 

646.7179 

0.0000 

1900 

3768.5610 

0.0000 

2200 

2628.0830 

0.0000 

2900 

437623.000 

0.0000 

3700 

1338404.00 

0.0000 
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Figure 3, Spectral absorption (top) and scattering (bottom) for oceanic constitu- 
ents in the oceanic component of OARM. Units for seawater absorption and 
scattering coefficeints are nr 1 , chlorophyll-specific absorption and scattering 
coefficients for phytoplankton groups have units of m 2 mg 1 . Note that seawater 
absorption coefficients have been scaled by a factor of 0.1, and seawater scatter- 
ing coefficients have been scaled by a factor of 10 to fit on the plot. 


Phytoplankton specific scattering coefficients are obtained from measurements at 590 
nm and extended to the entire spectrum from information on specific attenuation coeffi- 
cients (Bricaud et al., 1988; see Appendix). Again substantial group-dependent variabil- 
ity is apparent in the scattering spectra (Fig. 3). Additionally, knowledge of the specific 
backscattering coefficient is required for the oceanic model. These values are derived 
from the probability of backscatter (see Appendix). 

All calculations are performed in energy units (W nr 2 ), in conformance to the absorption 
and scattering derivations. When summing the mean irradiance in the oceanic layers of 
an Ocean General Circulation Model (OGCM), irradiance is converted to quantum flux 
(mmol quanta nr 2 s 4 ) to enable photoadaptation and growth kinetics to occur in accor- 
dance with the prevalence of laboratory analyses. 
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3 Discontinuities in the Cloud Data Set 


The ISCCP cloud cover and liquid water path data exhibited spatial discontinuities 
(Fig. 4). These discontinuities typically occurred in three general locations: the central 
Indian Ocean, western Atlantic Ocean, and the central Pacific Ocean. The Indian Ocean 
discontinuity was always the most obvious and persistent, the Pacific Ocean disconti- 
nuity was second, and the Atlantic Ocean discontinuity was diminished relative to the 
Indian and often not apparent. The discontinuities were always narrowest at the equa- 
tor and became wider toward high latitudes (Fig. 4). As such, they resembled the cov- 
erage typical of geostationary satellites. Similar observations have been noted by Norris 
(2000). These discontinuities were apparent in nearly every month and year of the data 
record, although the locations differed monthly on occasion and sometimes differed 
between the cloud cover and liquid water path data sets. 

The ISCCP data set is a compilation of radiances from many different satellites, but is 
generally derived from two series of geostationary and polar-orbiting missions. The 
radiances are combined and Level-3 products (geophysical variables such as the cloud 
cover and liquid water path of interest here) are derived and mapped to a standard Earth 
grid (Rossow et al., 1996). Although on-orbit failures and launch delays affect the 
coverage of the Earth by these two satellite series at any one time, there are usually 
three geostationary satellites, located in the western Pacific near 140°E (the GMS se- 
ries), in the eastern Atlantic near 0° (the METEOSAT series), and in the eastern Pacific 
near 90°W (the GOES series). Data from these three geostationary series are the prior- 
ity for the ISCCP-D2 archive (Rossow et al., 1996). Unfortunately, the GMS and 
METEOSAT satellites do not cover the central Indian Ocean, which is filled in by the 
NOAA polar-orbiting series. This geo-stationary gap lies over the Indian Ocean dis- 
continuity, immediately suggesting a cause. METEOSAT and GOES coverages over- 
lap in the central western Atlantic Ocean, corresponding to the second discontinuity. In 
cases of overlap, the ISCCP procedures are to preferentially use the data with the low- 
est solar zenith angle (Rossow et al., 1996). GOES and GMS overlap in the central 
Pacific again corresponding to the location of the Pacific discontinuity. 

Considering the effort devoted to inter-calibration among the diverse satellite sources 
by the ISCCP project (Rossow et al., 1996), we suspected that these discontinuities 
resulted from satellite/sensor inter-calibration inconsistencies that may have been too 
subtle to appear in the Level- 1 data but re-appeared in the Level-3 products as a result 
of extended processing. We attempted to ameliorate the effects of the sensor 
discontinuities in the Level-3 data, assuming that Level- 1 inter-calibration procedures 
by the ISCCP effort have been exhaustive. We first located the area of discontinuity in 
each monthly record. We eliminated all values that existed in the apparent discontinuities, 
by specifying a location in the center of the discontinuity at the equator, and its width, 
and then re-creating its shape in the mapped Level- 3 image by applying a cosine func- 
tion. The eliminated values produced gaps in the images of cloud cover and liquid 
water path (Fig. 4). 
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Figure 4. Depiction of the discontinuities in the ISCCP cloud cover (top 4 
panels) and liquid water path (bottom 4 panels) data sets for July 1989. This 
sequence shows the discontinuities in the upper left (note especially the cen- 
tral Indian Ocean), elimination of data occurring in these areas in the upper 
right, reconstruction using objective analysis in the lower left panel, and the 
final reconstruction where the edges of the reconstructed gaps are taken from 
objective analysis to serve as an interior boundary condition for the blended 
analysis. 









To reconstruct data in the newly created gaps, we first applied the Successive Correc- 
tion Method (SCM), using the Barnes method (Daley, 1991). We applied a radius of 
influence of 100 km. This method is often referred to as objective analysis. The SCM 
re-created data in the gaps in conformance with the spatial structure of the data outside 
the gaps produced by the geostationary satellites. However, it also produced a new 
discontinuity at the meridional center of the gaps where the spatial structure trends of 
the two geostationary satellites were in disagreement (Fig. 4). Thus we limited our 
application of SCM to a small number of gridded pixels inside the gap on each side. 
The number varied from 4-12 depending upon the size of the gaps, but 4 was used most 
often. Since these pixels are very close to the valid data provided by the geostationary 
satellites, we assumed these data to represent a “truth field”, and then applied the blended 
analysis of Reynolds (1988) to reconstruct the remainder of the gaps. This method 
produces a bias correction given a set of internal boundary conditions represented by 
the “truth field”, while adhering to the spatial structure of the original data in the gap 
area before elimination. This bias correction quality is the desired characteristic sought 
for this problem. The method has been used successfully in reconstructing sea surface 
temperature fields from satellite data using in situ observations (Reynolds, 1988, 
Reynolds el al., 1989), as well as for chlorophyll concentrations (Gregg and Conkright, 
2001). The resulting image fields indicated little apparent discontinuity, while retain- 
ing spatial structure of the original fields, even in the gap-eliminated regions (Fig. 4). 


4 Comparison of OAR M with Surface Solar Irradiance Data 

OARM results were compared to surface solar irradiance data produced by ISCCP 
(Bishop et al., 1997) to evaluate its performance. The ISCCP surface solar irradiance 
data represent monthly means integrated over the solar spectrum. OARM is a spectral 
and directional model over the same spectral domain. The comparison was performed 
only for the Jan. 1989-May 1991 using coincident monthly mean atmospheric optical 
variables (not the climatologies discussed earlier). This time period was selected be- 
cause of the coincident availability of ISCCP cloud data from the GES-DAAC and 
TOMS ozone data. Although the ISCCP surface solar irradiance data represented the 
output from a model (Frouin et al., 1989), it has been tested extensively against field 
measurements with excellent results (Bishop et al., 1997). The comparison performed 
here involved integrating the surface spectral irradiance computations of OARM across 
the entire solar spectrum and over the direct and diffuse components, with surface re- 
flectance neglected. Thus we effected a one-to-one comparison with data, that indi- 
cated the performance of the spectral and directional OARM. 

The comparison between OARM integrated spectral and directional irradiance exhib- 
ited a large degree of correspondence with ISCCP surface solar irradiance for the 29 
month record under investigation (Table 4). Mean coefficient of determination (r 2 ) was 
0.97, with RMS of 12.1% and an average error of -4.0%. Minimum and maximum r 2 
was 0.96 and 0.98, respectively. The minimum RMS was 9.9% occurring in Jun. 1989 
and the maximum was 15.1% occurring in Mar. 1990. The maximum was primarily the 


11 



result of disagreements occurring in the Antarctic Ocean. The minimum average error 
was -7.3% and the maximum was -1.2%. The negative sign on the average error indi- 
cates bias: OARM was typically lower than ISCCP. 

Table 4. Comparison of OARM surface solar irradiance integrated over the solar 
spectrum and over direct and diffuse components with total surface solar irradiance 
data from ISCCP over a 29-month period from Ian. 1989 to May 1991. 



Mean 

Minimum 

Maximum 

r 2 

0.97 

0.96 

0.98 

RMS 

12.1% 

9.9% 

15.1% 

Average Error 

-4.0% 

- 7.3% 

-1.2% 


Representative scatterplots of OARM and ISCCP total surface irradiance for four sea- 
sons indicated the correspondence between the two (Fig. 5). The vast majority of the 
points in each month fall on a one-to-one line, but OARM had a substantial number of 
lower computed irradiances in the mid-range (about 120-250 W nr 2 ) than ISCCP. Out- 
liers well below the one-to-one line indicated differences in ice masking data sets used 
here: ISCCP used a radiance threshold (Bishop et ai, 1997) while OARM used ice 
sheet distributions from the NASA Data Assimilation Office. In Feb. 1990, Aug. 1990, 
and Nov. 1989 the comparison indicated slightly higher estimates by OARM at low 
irradiance levels (< 50 W nr 2 ) than ISCCP. This was due to Earth curvature effects that 
were accounted for in OARM, but not in ISCCP. 
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Distributions of integrated surface irradiance for Feb. and Sep. 1990 (Figs. 6 and 7) 
i llustrated the nature of the agreement between OARM and the ISCCP data sets. Gen- 
eral features of surface irradiance were in agreement throughout the global domain, 
especially over the oceanic portions. The ISCCP data exhibited smoothing of the re- 
trieved radiance fields (J.K.B. Bishop, personal communication, 1999). The input cloud 
fields used by OARM were not smoothed. This gave rise to some of the differences 
observed in the statistical analyses. Overall, these distributions of surface integrated 
irradiance reflected the correspondence observed in the statistical analyses. 
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Appendix 

Mathematical description of the oceanic component of OARM 


The oceanic portion of OARM tracks three irradiance streams in the water column: 
downwelling direct, downwelling diffuse, and up welling diffuse. We track only one 
scattering event for each of the oceanic layers for computational convenience. 

The direct downwelling stream is modeled as 

E d (X,k) = E d (/.,kri)T d (A.k)T d (Ak) (Al) 

where E d is the downwelling direct spectral irradiance at level k in the water column 
(corresponding to the layer definitions of an OGCM, where k = 1 is the surface layer, 
and k is positive downward), T ad is the spectral transmittance of the direct stream after 
absorption, and T gd is the transmittance after scattering 


T d = exp(-a/p d (0 w ) Az) 

(A2) 

T^ = exp(-b/p d (0 w ) Az) 

(A3) 


(^-dependences have been dropped), where a is the total absorption coefficient of the 
water in OGCM layer k, b is the total scattering coefficient, u d (0 w ) is the average co- 
sine, or path length in the water column as a function of the downwelling zenith angle 
0 w , and Az is the thickness of the OGCM layer in question. The downwelling zenith 
angle is related to the solar zenith angle 0 by the index of refraction of seawater 


U.(0 ) = COS0 

d v w 7 w 

(A4) 

= sin 4 [sin0/n] 

(A5) 


where n is the seawater refractive index, parameterized as 1,341 (Austin, 1974). 

The downwelling diffuse irradiance gains forward scattered downwelling direct irradi- 
ance and backscattered upwelling diffuse irradiance in the water column 

E s (k) = E..(k-1 )T as (k)T ss (k)+F[ I -T ss (k) j + E d f k-l)T d (k){F[ 1 -TJk)]} 

+ E u (k+ 1 )T au (k) { B [ 1 -T su (k)] } (A6) 

where E is the downwelling diffuse irradiance at OGCM laver k, T and T are the 
diffuse transmitt ances after absorption and scattering, respectively, F and B are the 
probabilities of forward and back scattering, respecti vely, T au and T su are the upwelling 
diffuse transmittances after absorption and scattering, respectively, and E u is the up- 
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welling diffuse irradiance. T ,, T,., T , and i are defined similarly to their down welling 
direct counterparts, except that p d is replaced by p s for downwelling direct irradiance, 
and by p u for upwelling diffuse irradiance. They are considered constant in OARM at 
0.83 (Sathyendranath and Platt, 1988) and 0.4 (Kirk, 1981), respectively. The first 
term on the right-hand side of Eq. A5 represents transmittance after scattering and ab- 
sorption for the downwelling diffuse stream, plus the addition of the forward-scattered 
component. The second term represents forward scattering from the downwelling di- 
rect stream from the layer above, and the third term represents backscattered irradiance 
from the upwelling diffuse stream from the layer below. The backscattering probabil- 
ity B is defined as 


B = b b /b (A7) 

where b b is the backscattering coefficient. The forward scattering probability is simply 

F = 1 - B (A8) 

The diffuse upwelling stream is derived from the two-stream model of Aas (1987), with 
contributions from both the downwelling direct and downwelling diffuse explicitly in- 
cluded, similarly to Ackleson et al. (1994). After experimentation, we found that the 
average cosine for downwelling diffuse irradiance could be substituted for the direct 
average cosine without incurring > 1 pmoles quanta nr 2 s" 1 error, except in the compu- 
tation of irradiance leaving the ocean surface . With these assumptions, the integral form 
of the Aas (1987) model reduces to 


: a (k) = [E d (k-1)+E s (k-1)] exp(cxAz) (a+c s )/b a 

(A9) 

0.5 {(c u -c s ) - [(c s -c u ) 2 -4(b s b a -c s c u )] 1/2 } 

(A10) 

c s = (a+rb .)/ p s 

(All) 

c a = (a+r a b b )/ fi a 

(A12) 

K = r s b h /u s 

(A13) 

K = r A/p» 

(A14) 


where r s and r u are dimensionless shape factors to represent the mean backscatter coef- 
ficient (Aas, 1987), and are assumed constant at 1.5 and 3.0, respectively (Ackleson et 
al., 1994). 

Phytoplankton specific absorption coefficients a* p (X) are derived by taking spectra from 
the published literature and normalizing to the absorption at 440 nm [a* p (440)] to ob- 
tain normali zed specific absorp ti on spectra [a* (X)] N for each of the five phy toplankton 
groups. Diatom and chlorophyte [a* p (7.)] N are taken from Sathyendranath et al. (1987), 
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cyanobacteria from Bricaud et al. (1988), coccolithophores from Morel and Bricaud 
(1981), and dinoflagellates from Ahn et al. (1992). Then the specific spectral a* p (>.) 
values are derived using mean values at 440 nm. Diatom a* p (440) represents the mean 
of 5 observations containing 4 different spp., chlorophytes 6 observations from 4 spp., 
cyanobacteria 5 observations from 3 spp., coccolithophores 3 observations of 1 spp., 
and dinoflagellates 1 spp, all from the references listed above. CDOM absorption is 
taken from Carder et al. (1989) 

a CDOM (X) = c h a* h (X) + C f a * { (k) (Al 5) 

where a CDOM is the absorption coefficient of CDOM (m 4 ), C represents the concentra- 
tion (mg nr 3 ) of a humic acid component (subscript h) and a fulvic acid component 
(subscript f), and a* is the specific absorption coefficient of the two components 

■d* h (k) = 1 .302X10- 4 exp[450-(X)] (Al 6) 


and 


a* f (l) = 5 .0x1 O' 6 exp[450-(X)] (A17) 

(Carder et al., 1989). C h and C f are fixed as 1 .0 and 10.0 mg nr 3 , respectively (Ackleson 
et al., 1994). 

Phytoplankton specific scattering coefficients are obtained from measurements at 590 
nm and extended to the entire spectrum from information on specific attenuation coef- 
ficients (Bricaud et al., 1988). Diatom and clilorophyte specific scattering coefficient 
at 590 nm, b* d (590) and b* c (590), are the mean of 5 observations and 6 observations, 
respectively, from Morel (1987), Bricaud and Morel (1986) and Bricaud et al. (1988). 
Cyanobacteria b*.(590) is the mean of 8 observations from Morel (1987), Bricaud and 
Morel (1986), Bricaud et al. (1988) and Ahn et al. (1992). Coccolithophore b* o (590) is 
derived from the mean of 3 observations from Bricaud and Morel (1986), Bricaud et al. 
(1988), and Ahn et al. (1992). Dmoflagellate b* { (590) is derived from a single observa- 
tion by Ahn et al., (1992). 

Measured chlorophyll-specific attenuation coefficients (Bricaud et al., 1988) suggested 
1) a slight negative linear spectral trend for diatoms, 2) a spectrally flat profile for 
chlorophytes and coccolithophores, 3) a slight negative trend for cyanobacteria, with a 
local maximum corresponding to phycoerythrin absorption. But Bricaud et al. (1983) 
observed a strong spectral dependence for Emiliana huxleyi, which we adopt here. Given 
this information and knowledge of c* p (590) for each phytoplankton group 

c* p (590) = a* p (590) + b* p (590) (Al 8) 
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we fit the total attenuation spectra of Bricaud et al. (1988) by 


c* d (X) = c* d (590) - 4.3xl0- 4 (X-590) 

(A19) 

c* c (X) = c* c (590) 

(A20) 

c * (X) = c* (590) - 8.0xlO 4 (X-590) 

(A21) 

c* 0 (X) = c* o (590) - 8.95xl0 4 (X-590) 

(A22) 

c* f (X) = c* t (590) 

(A23) 


assuming no spectral dependence for dinoflagellates. 

The specific total backscattering coefficients are derived by subtraction 

b*(X) = c*(X) - a*j(X) (A24) 

where the subscript i indicates the phytoplankton group of interest, diatoms (d), 
chlorophytes (c), cyanobacteria (y), coccolithophores (o), and dinoflagellates (f). 


Additionally, knowledge of the specific backscattering coefficient is required for the 
oceanic model. These values are derived from the probability of backscatter 

b b (X) = B x b(X) (A25) 

We assume no spectral dependence in b b> and use a constant value of 0.002 for diatoms 
(Morel, 1988), 0.00071 for chlorophytes, 0.0032 for cyanobacteria (Ahn et al. 1992), 
0.00071 for coccolithophores, and 0.0029 for dinoflagellates (both from Morel, 1988). 
Ahn et al. (1992) suggested a spectral dependence for cyanobacteria but generally none 
for the other groups. 

Phytoplankton respond to the mean irradiance availability in the isopycnal layers in 
which they reside. This is especially important in an OGCM simulation with discrete 
numerical layers. OARM computes irradiance at the vertical boundaries of the model 
layers. Derivation of accurate mean irradiance within each model layer is essential for 
simulating realistic phytoplankton responses and primary production estimates in the 
context of variable cloud cover and thickness. 


We approach this problem by deriving a lookup table of factors to convert irradiance at 
boundaries to the mean within the layer. The lookup table incorporates depth, chloro- 
phyll concentration, and phytoplankton group as array indices. CDOM is assumed 
constant, but variability could be incorporated by addition of another dimension in the 
lookup table. The table is created applying the oceanic model in advance with the 
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assumption that p d = 1 uniformly. Computation of mean irradiance is 

E t avg(z) = M(Az,chl,p) { [E t (z- 1 )+E t (z)]/2 } (A26) 

where M is a factor obtained from the lookup table incorporating layer thickness (Az), 
chlorophyll concentration (chi), and phytoplankton group (p), and E t is the total irradi- 
ance. Note that this computation is non-spectral, although the computation of irradi- 
ance at the layer boundaries is spectral. The effect of the computation for average 
irradiance is not subtle - lookup table factors range from 0.0009 at high chlorophyll 
concentrations and 200 in layers to 0.95 at low chlorophyll and shallow layers (< 5 m). 
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